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Abstract. The atomic cascades in /i~p and pp atoms have been studied in detail using new results for 
the cross-sections of the scattering of highly excited exotic atoms from molecular hydrogen. The cascade 
calculations have been done with an updated version of the extended standard cascade model that computes 
the evolution in the kinetic energy from the beginning of the cascade. The resulting X-ray yields, kinetic 
energy distributions, and cascade times are compared with the experimental data. 

PACS. 36.10.-k Exotic atoms and molecules (containing mesons, muons, and other unusual particles) 



1 Introduction 

The standard cascade model (SCM) of exotic hydrogen 
atoms, originally introduced by Leon and Bethc [Q and 
later extended to include the evolution of the kinetic en- 
ergy distribution during the atomic cascade [^,|^] , provides 
a fair description of many properties of atomic cascades, 
such as the X-ray yields and absorption fractions M fj\. 
However, no detailed calculations of the initial stage of 
the cascade were done until recently due to the lack of 
the corresponding cross-sections. The upper stage of the 
cascade was commonly described with a phenomenologi- 
cal deexcitation mechanism, the so-called chemical deex- 
citation This mechanism was introduced in order to 
explain the data on cascade times, as the external Auger 
effect was found to be too slow at the initial stage of the 
cascade. The exact nature of the collisional deexcitation of 
highly excited exotic atoms remained a mystery for a long 
time until experimental studies of the initial stages of the 
atomic cascade became possible. Recently the energy dis- 
tributions of n~p, fJ-~p, and fi~d atoms were studied with 
various time-of-flight methods P~p^|. In particular, the 
measurements performed with fi~p and pp at low density 
allow one to probe collisional deexcitation mechanisms for 
highly excited states as long as the energy distribution re- 
mains frozen during the lower cascade stage dominated by 
the radiative transitions. 

In our recent paper [DJ, we studied the dynamics of 
collisional deexcitation of highly excited exotic atoms us- 
ing the classical-trajectory Monte Carlo method. The Coulomfthe exotic atom: the kinetic energy distribution changes 
transitions with large change of principal quantum num- during the cascade due to the acceleration and deceler- 
ber n were found to be the dominant collisional deexcita- ation mechanisms ^, [l4|, [l5|] . The new results for the 
tion mechanism at high n, with the molecular structure of collisional processes presented in Jl3|,[l^] are used in the 
the hydrogen target being essential for the dominance of present version of the ESCM, a significant improvement 
transitions with large An. The main goal of this paper is over the previous calculations. 



to investigate the atomic cascades in hydrogen-like atoms 
at low density and to confront the theoretical results with 
the experimental data on the X-ray yields, cascade times, 
and kinetic energy distributions for muonic and antipro- 
tonic hydrogen. As the lower part of the atomic cascade 
at low target densities is mainly dominated by the radia- 
tive deexcitation, the results of our calculations are only 
weakly affected by theoretical uncertainties in the colli- 
sional cross-sections for the low n states. The detailed cas- 
cade calculations for high density targets will be published 
elsewhere. 

The paper is organized as follows. The cascade model is 
described in Section 0. The results of the cascade calcula- 
tions for the u~p &napp atoms are presented in Section^. 
The conclusions are summarized in Section 0. 



2 The extended standard cascade model 
2.1 Overview 

The extended standard cascade model (ESCM) is a ki- 
netics model that includes all cascade processes of the 
standard cascade model Jl|,§,0: the radiative, Auger, and 
Coulomb deexcitation, Stark mixing, and, in case of hadronic 
atoms, nuclear absorption (see Tableland references therein). 
In addition, the ESCM also takes into account the inter- 
action between internal and external degrees of freedom of 
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Stark mixing: 
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External Auger effect: 
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H -> 
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Coulomb deexcitation: 
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~P)ni 


f Ha 


-> (z~p) n/ +H2, n f < n,i 
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Elastic: 
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f H 2 


-> (x~p)„i + Hj 
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Absorption: 
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Radiative: 
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~p)n ( l, 






Nuclear reaction: 


(7T 


~p)ns 


-> V 


+ n, 7 + n 


[i] 


Weak decay: 


7T 


-> /x" 









Table 1. Processes included in the extended standard cascade model. 



The cascade in exotic atoms is divided into two parts: 
the classical domain for high n and the quantum mechani- 
cal domain for low n. The x~p can be described classically 
for quantum numbers n>l, and for convenience we use 
the properties of the Auger deexcitation to define the n- 
ranges of the two domains. The Auger deexcitation rates 
are known to have a maximum for n = n c where n c is the 
largest n for which the binding energy released in a An = 1 
transition is enough to ionize the H 2 molecule @,[[j|Jl|: 
n c = 7 for n~p and n c = 12 for pp. Our calculations show 
that the Coulomb deexcitation dominates the Auger de- 
excitation for n > n c |l3| . We define the classical domain 
by the conditions n > n c . The classical domain, where 
the processes were calculated in the classical Monte Carlo 
(CMC) model with the molecular structure of the target 
taken into account, will be discussed in Section 2.2. In the 
quantum mechanical domain, n < n c , the close-coupling 
method and semiclassical approximations Jl6|,[ij],|2l[] were 
used to calculate the differential cross-sections, der/d cos 8, 
of the collisional processes 



(x p)nl + H -> (x p) n l> + H 



(1) 



and the cross-sections for absorption during collision and 
Auger transitions 



(x p) n i + H 



absorption, 
{x~p)n, +P 



(2) 
(3) 



in an energy range relevant for atomic cascade. The new 
calculations allow us to describe the competition between 
deceleration and Stark mixing, as well as the absorption in 
hadronic atoms due to Stark collisions, without employing 
any fitting parameters related to the Stark mixing and 
deceleration, like fcgtk used in many earlier calculations 0). 
Our cascade code makes full use of the differential cross- 
sections so that the kinetics is treated more accurately 
than in cascade models using continuous deceleration |g2| . 

The differential cross sections used in the present cas- 
cade model in the case of muonic hydrogen were calculated 
in the fully quantum mechanical close-coupling framework 
for n = 2 — 5 and in the semiclassical approximation for 
n = 6 — 7. The statistically weighted differential cross- 
sections were used for n > 6. In antiprotonic hydrogen, 
the fully quantum mechanical results are not yet available 
below the ns thresholds. Therefore, we used the results of 
the semiclassical model in the range n = 2 — 9 and the 
fixed field model for n = 10 — 12. 



Another approximation used in this paper is related 
to the Auger deexcitation because the eikonal approxima- 
tion |l^] does not give the differential cross-section and 
the distribution over final I. Here we use the differential 
cross-section for the Z-average Stark and elastic transitions 
and a statistical distribution over final I. 

Concerning the acceleration mechanisms, the present 
calculations include Coulomb deexcitation through the whole 
cascade. For high n states, the classical-trajectory Monte 
Carlo method |L3|] was used to obtain the cross-sections of 
the inelastic collisions 

(x-p) ni + H 2 - (x~p) nf + X, X = H 2 , H*, H + H.(4) 

For low n, where the classical-trajectory method cannot 
be expected to give reliable results, we parameterize the 
Coulomb deexcitation cross-section as follows[] 



i-i(E c 



(5) 



where E cm is the CMS kinetic energy and and fi xp are 
the reduced masses of the pTp and the x~p, respectively. 
We use the values 



c = 1.2 • 10' 3 al eV, 
7 = 3.5 



(G) 
(7) 



which gives a fair description of the Coulomb cross-sections 
of Bracci and Fiorentini [23j] for n < 7. 



Other calculations J24 - |26| predict significantly smaller 
Coulomb deexcitation rates than |23| . However, the rapid 
increase in the Coulomb cross-sections for energies ap- 
proaching zero is predicted by all models. The dependence 
on the reduced mass is based on an estimate in [Q and 
our classical Monte Carlo calculations confirm that it is 
a fair approximation. We assume, furthermore, that only 
An = 1 transitions are important at low n, the distribu- 
tion over final / is statistical, and the angular distribu- 
tion is isotropic. As long as the mechanism responsible for 



The n = 2 state in muonic hydrogen is, however, treated 
differently. Coulomb deexcitation is assumed only to take place 
in the 2s state and only when the kinetic energy is below the 
2p threshold. In this case we use the measured Coulomb de- 
excitation rate Ag^ff b = 4.4 • 10 n iV s^ 1 ft| where N is the 
density of the target in units of liquid hydrogen density (LHD). 
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Fig. 1. The /-average rates at T = 1 eV for muonic hy- 
drogen in gaseous target at 10 -4 LHD. The Coulomb de- 
excitation (filled diamonds) and Stark mixing (filled trian- 
gles) rates calculated in the classical-trajectory model (l]] are 
shown in comparison with the results of the semiclassical fixed 
field model for Stark mixing (solid line), Auger deexcitation 
(light dashed line), and nuclear absorption during collisions 
(dashed line) |f^,Q. The Coulomb deexcitation rate from [^3| 
is shown with a dash-dotted line. The radiative np — » Is and 
n(n — 1) — » (n — l)(n — 2) rates are shown with dotted lines. 



Fig. 2. The Z-average rates at T = 1 eV for antiprotonic 
hydrogen in gaseous target at 10~ 4 LHD. The Coulomb de- 
excitation (filled diamonds) and Stark mixing (filled triangles) 
rates calculated in the classical-trajectory model [[[3) are shown 
in comparison with the results of the semiclassical fixed field 
model for Stark mixing (solid line), Auger deexcitation (light 
dashed line), and nuclear absorption during collisions (dashed 
line) [ fn^[l6[ . The Coulomb deexcitation rate from Equation ^ 
is shown with a dash-dotted line. The radiative np — > Is and 
n(n — 1) — > (n — l)(n — 2) rates are shown with dotted lines. 



Coulomb deexcitation at low n is not fully understood^, 
this process enters as a major uncertainty in calculations 
of the kinetic energy evolution at low n. In this paper we 
restrict the cascade calculations to observables that are 
not very sensitive to the Coulomb deexcitation at low n. In 
the case of kinetic energy distributions this usually means 
low densities. Other observables, like the X-ray yields in 
H~p and cascade times, are less sensitive to the Coulomb 
deexcitation at low n. 

Figures [l] and || show an overview of the cascade rates 
at 10~ 4 LHD in muonic and antiprotonic hydrogen calcu- 
lated in the classical-trajectory model with molecular tar- 
get and in semiclassical approximations with atomic tar- 
get. The absorption rate in antiprotonic hydrogen shows 
only the absorption from the ns states during collisions. 
Cascade calculations show, in agreement with the exper- 
imental results p9[ , that absorption at densities below 
10~ 3 LHD takes place mainly from the p states. 



2.2 Classical-trajectory Monte Carlo model in the 
cascade calculations 

The classical-trajectory Monte Carlo calculations have been 
included in the cascade calculations by using two different 
methods. The earlier versions of the cascade model [14| 
was based on the explicit use of the cross-sections [13| 
with the following initial conditions. A discrete set of 9 
initial kinetic energies in the interval 0.05 eV < T < 20 eV 
was used, and the initial principal quantum numbers were 



taken in the range 8 < n.; < 20 for [i~p and 13 < n.; < 35 
for pp. For each combination of laboratory kinetic energies 
T and principal quantum numbers rii, the classical trajec- 
tories were calculated using the statistical distribution in 
quantum number li as the initial condition. The differ- 
ential cross-sections dcr/dcos6> were calculated separately 
for the three groups of reactions: 



riili elastic 
n,*// (If ^ li) Stark 
riflf (rif < rii) Coulomb 



(8) 



2 The investigation of the role of molecular resonances [ p7|[28[ 
is beyond the scope of this study. 



The effect of excitation and dissociation of the hydrogen 
molecule in Coulomb deexcitations were taken into ac- 
count by distributing the final kinetic energies of the three 
atoms (x~p and two H) in the CMS according to the clas- 
sical phase space. This approximation agrees fairly well 
with the calculations [jl3| . 

More detailed initial and final states are used in the 
present cascade model where several trajectories (21 for 
\x~p and 6 for pp) were calculated for each combination 
of the quantum numbers rii (in the range given above) 
and h — 0, — 1, and 31 kinetic energy points in the 
range 0.3 — 32 eV. In muonic hydrogen this gives a total 
of 118482 trajectories. The range of impact parameters 
was divided into three intervals, (0;2ao), (2ao;4ao), and 
(4ao;7eto), with a third of the collisions in each and the 
trajectories uniformly distributed in p 2 . For each trajec- 
tory, the final quantum numbers nf and If, and the final 
laboratory kinetic energy Tf were stored. The resulting 
lists of transitions in the form 



Til li Ti 



(9) 
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with the weights of the impact parameter intervals prop- 
erly taken into account, were used directly by the cascade 
program by redirecting the final state to the initial state of 
the next collision. We consider this method as more accu- 
rate than the earlier used method for cascade calculations 
in the upper part of the cascade. 



2.3 Initial conditions 

The initial conditions for the cascade calculations are de- 
fined by the initial distributions in the quantum numbers 
n and I and the laboratory kinetic energy T of the x~p. 
In the simplest picture, the x~ is captured by the pro- 
ton after ionization of the hydrogen atom in a state with 
large overlap with the Is electron wave function; this gives 
Tii ~ y Hxp/me (~ 14 for muonic and ~ 30 for antipro- 
tonic hydrogen). More elaborate calculations, which takes 
the molecular structure of the target into account, predict 
distributions in the initial n which peak at lower values: 
for example, the maximum at n = 11 for muonic hydro- 
ge n M . A shift towards lower initial n values is also found 
in [|3 lj] 7 The new experimental data on muonic hydrogen 
kinetic energy distributions from PSI are sensitive to the 



initial distributions as discussed below in Section 3.1.4 



Unless otherwise indicated, we use the initial principal 
quantum number ni = 14 for fi~p and rii = 30 for pp and 
the initial kinetic energy given by a Maxwell distribution 



At target densities above 10 -7 LHD, the lifetime of the 
2s state depends on the kinetic energy of the ji~ p. Below 
the 2p threshold, T 2p = 0.3 eV, the 2s state is long-lived 
because 2s — > 2p Stark transitions followed by the fast 
radiative deexcitation 2p — > Is are energetically forbid- 
den. The induced radiative transition can, however, occur 
during collisions, see jy],|^] and below. For kinetic ener- 
gies above the 2p threshold the collisions with the target 
molecules lead to a competition between deceleration and 
depletion via 2s — * 2p Stark transitions followed by the 
2p — > Is radiative transition. The metastable 2s fraction, 
i?2s, is defined as the fraction of all created [i~p atoms 
which end up in the 2s state with kinetic energies below 
T 2p . 

Before we present the results of our full cascade model, 
it is illuminating to consider the problem of the metastable 
2s fraction in a simplified framework. The metastable 2s 
fraction consists of the fi~p that are formed with kinetic 
energies below T 2p and those that are formed with kinetic 
energies above T 2p but slow down below threshold due to 
elastic collisions. The median kinetic energy at 1 mbar is 
1.5 eV Q, so the deceleration process is important. This 
leads to the definition of the surviving metastable fraction, 
f{T), which is the probability that the fi~p atom in the 2s 
state with initial kinetic energy T reaches energies below 
the 2s threshold by slowing down in elastic collisions |2lJ . 

In the approximation of continuous energy loss the evo- 
lution in the kinetic energy is given by 



w(Ti) — 3\j J^, ex P 



2ttT 



2ttT 



(10) 



with To = 0.5 eV. In all cases the k distribution is statis- 
tical. 



3 Results 



T — 

dt 



Mb 



rNva% s _ 



-AT) (11) 



where M pp and Mh are the masses of the fi~p and the H, 
respectively, v the velocity of the fi~p, and 02s->2sO^) the 
transport cross-section 
Carboni and Fiorentini 
of the surviving metastable fraction 



. Equation (Jill) was used by 
to get the following estimate 



3.1 Muonic hydrogen 

The cascade calculations in muonic hydrogen have been 
done in the density range (10 -8 — 1) LHD using various 
initial conditions. The results are compared with the ex- 
perimental data for the metastable 2s population, X-ray 
fields, cascade times, and kinetic energy distributions jl2], 
101. 



3.1.1 The metastable 2s fraction 

The metastable 2s state of muonic hydrogen is of great in- 
terest because it allows, under certain conditions, to mea- 
sure the \i~p Lamb shift, 2s — 2p, from which the proton 
charge radius can be determined with high precision Q . 
The feasibility of this experiment, presently in progress at 
PSI p9[|, relics on a sizeable fraction of muonic hydrogen 
atoms in the 2s state with a sufficiently long lifetime. In 
the absence of collisions, the lifetime of the 2s state is de- 
termined by the muon lifetime (t m = 2 • 10~ 6 s) as the rate 
for the 2s — > Is two-photon transition is negligibly small. 



f(T) = exp 



(M w + Af H ) s 
2M W M H 



AT 



T 2p 1 a 2s- 



-dT' 



(12) 

It was found that a sizeable fraction of (/z p) 2s atoms 
formed at kinetic energies below 1.3 eV slows down below 
threshold. 

To provide a more realistic treatment of the evolution 
in kinetic energy we use a Monte Carlo program based 
on the differential cross-sections for the four processes 
2s — + 2s, 2s — > 2p, 2p 2s and 2p — > 2p. In addition 
to the collisional processes, the 2p — > Is radiative transi- 
tion is also included in the calculation. The fate of a [i~p 
formed in the 2s state with kinetic energy T is thus either 
to undergo 2p — + Is radiative transition after the Stark 
mixing 2s — > 2p or to end up in the 2s state with kinetic 
energy below the threshold with probability f(T). 

Figure U shows the surviving metastable fraction, f(T), 
calculated with the Monte Carlo program for target den- 
sity 10" 6 LHD < N < 10" 2 LHD. The approximation (jl|) 
gives somewhat higher values for the survival probability 
than the exact kinetics calculation at T < 1.4 eV. The 



T.S. Jensen, V.E. Markushin: Collisional deexcitation of exotic hydrogen atoms. II. Cascade calculations 



5 



Monte Carlo results at high energies (T > 1.5 eV) are sig- 
nificantly larger than those obtained from Equation ( |l2| ) 
where continuous energy loss is assumed. The reason is 
that there is a significant contribution to f(T) from large- 
angle scattering. The result of |h| is considerably smaller 
than the present result due to differences in the cross- 
sections and the fact that in (io| the molecular hydrogen 
mass was used instead of the atomic one in order to get 
the correct value for the threshold. By using the surviving 
metastable fraction /(T), the measured kinetic energy dis- 
tributions on arrival in the Is state [ p2| , and the 2s arrival 
probability^ [Q, one finds the metastable 2s fraction to 
be ~ 1% in the pressure region 0.06 — 16 mbar as discussed 
in & 

After reaching thermal energies the (// p) 2s atoms may 
undergo radiative deexcitation during collisions because of 
the Stark effect 

(n-p)2s+H -» {fi-p) 2 s-2 P +ii -» (/x"p)i s +7(1.9 keV)+H. 

(13) 

The radiative quenching of thermal (/i p) 2s atoms would 
result in delayed K a X-rays. However, the experiments |32| - 
HH have not observed this effect. This apparent contra- 
diction was recently resolved by the \i~p diffusion experi- 
ment [ [l2| , |43| ] where 0.9 keV (/i~p)i s atoms, resulting from 
the Comomb deexcitation process 

(ji-p) 2 . + H -> (m~p)i.(0.9 keV) + H(1.0 keV), (14) 

were found. The measured non-radiative quenching rate 
is about 2 /xs _1 at 1 mbar |l2|,^| which is more than an 
order of magnitude larger than the theoretical predictions 
for the radiative quenching rate j4l], [l2| . 

The 2s arrival probability, e 2s , and metastable 2s frac- 
tion calculated in th e E SCM with the initial conditions 
specified in Section 2J3 is shown in Figure ||. The cal- 
culated 2s arrival probability is in good agreement with 
the experimental data obtained from measurements of X- 
ray yields p2]-|34|]. The arrival probability e 2s increases 
steadily with increasing density from 1% at 10~ 7 LHD 
to 25% in liquid hydrogen. This behavior has a simple 
qualitative explanation: at the lowest densities the cas- 
cade proceeds purely radiatively from high n values. As 
only Al = 1 transitions are allowed, the circular states 
are predominantly populated for low n. The majority of 
exotic atoms, therefore, go through the radiative 2p — > Is 
transition without reaching the 2s state. With the density 
increasing the Stark mixing becomes efficient at lower n 
levels and more np — > 2s radiative transitions take place. 
In liquid hydrogen, e 2s is given by the statistical weight of 
the 2s state. 

In the density range 10~ 7 — 10~ 4 LHD, the calcu- 
lated metastable fraction increases from 0.7% to 2%, and 
this is in perfect agreement with the experimental re- 
sult [JfJI obtained in measurements of 0.9 keV ([i~p)i s 



1.0 



Detailed kinetics 



3 The 2s arrival probability is the fraction of formed exotic 
atoms that arrives at the n — 2 level in the 2s state. It can 
be estimated from measurements of K X-ray yields using the 
known Y„ p ^2s /Y np ->i s branching ratios. 
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Fig. 3. Energy dependence of the surviving metastable frac- 
tion f(T) of the 2s state of muonic hydrogen. The result of the 
detailed kinetics calculations (solid line) is shown in compar- 
ison with the approximation of continuous deceleration ( p^| 
(dashed line) and the result of Carboni and Fiorentini E| . 
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Fig. 4. Ab initio cascade calculation of the 2s arrival prob- 
ability, €2s (dashed line), in comparison with the values de- 
rived from X-ray measurements: (circles), ^] (triangles), 
and [^3) (filled square). The solid line shows the the metastable 
2s fraction, R,2 S ; the corresponding experimental data (filled 
diamonds) are from [|L2|. 



atoms. The metastable fraction grows slowly with den- 
sity until about 0.01 LHD. Above 0.01 LHD the growth 
becomes faster and R 2s reaches 65% in liquid hydrogen. 
The large metastable fraction at high densities is a conse- 
quence of the dominance of the slowing down in 2U — > 21 f 
collisions over the radiative 2p — > Is deexcitation. 

Figure || shows the time dependence of the 2s popu- 
lation, P 2s , at 1 mbar. The 2s population grows rapidly 
with time in the beginning of the cascade and reaches 
a maximum of 1.2% after 0.1 /is. The development at 
later times clearly shows the two components of P 2s . The 
short-lived component consists of the ^Tp atoms with ki- 
netic energies above the threshold, T 2p , which go through 
the 2s — > 2p Stark transition followed by the radiative 
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Fig. 5. The time dependence of the 2s population, P23 
(solid line), of muonic hydrogen at 1 mbar. The thin dashed 
line shows P2 S calculated without collisional quenching of the 
metastable 2s state. The dashed-dotted line shows the distri- 
bution of the K a X-rays in units of /is" 1 . 



deexcitation. The long-lived component consists of the 
metastable (u~p) 2s - The atomic cascade is essentially ter- 
minated after 0.7 /is: the u~p atoms are either in the Is 
state or the metastable 2s state. The metastable 2s state 
is quenched non-radiatively ( |l4|) with a rate at 1 mbar of 
about 4 times the rate for muon decay. 

The results in Figure || were calculated under the as- 
sumption that all \jTp atoms were formed at the same 
time. In reality, the muons are injected into the target, 
slow down in collisions with the H2 molecules, and are then 
captured by a proton. Having to go through the slowing 
down process means that the time of formation is given 
by a density dependent distribution — the so-called stop 
time distribution. The measured stop time at 1 mbar is 
Tstop = 0.48±0.13 fis |l2| which, when taken into account, 
does not change the overall features of Figure || except for 
the time offset. 

The results for the 2s population are extremely impor- 
tant for the 2s — 2p Lamb shift experiment p^,p9fl. In this 
experiment, a laser will be used to induce 2s — > 2p tran- 
sitions which arc followed by the fast radiative 2p — > Is 
deexcitation. The K a X-rays are detected in coincidence 
with the laser pulse and the 2s — 2p energy splitting is 
then deduced from the laser wavelength corresponding to 
the strongest K a signal. This laser experiment is feasible 
provided the 2s population is not too small and the K a 
background produced by the cascade is negligible. These 
conditions are fulfilled at very low densities (around a few 
mbar) after a short delay time as discussed above. 



3.1.2 X-ray yields 

The calculated relative X-ray yields for the fi~p atom are 
shown in comparison with the experimental data in Fig- 
ure ||. The agreement between theory and experiment is 
good for the K a , Kp, and yields. The relative K a 
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Fig. 6. The density dependence of the relative X-ray yields, 
K a , Kp, and K 1 , and the absolute total yield Ktat in muonic 
hydrogen. The experimental data are from 
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Fig. 7. The ratio of the K a and Kp yields for the non-radiative 
(solid line) and radiative (dashed line) quenching of the 2s 
state. The data are from ■ 



yield is higher than 90% at 10 7 LHD, it decreases with 
the density increasing and reaches a minimum of 50% at 
2 • 10~ 3 LHD. Above 2 • 10" 3 LHD the relative K a yield 
increases and ends at 94% in liquid hydrogen. The op- 
posite behavior is seen in the other K yields. They start 
out close to zero at 10~ 7 LHD, increase with the density 
and reach their maxima (at 0.01 LHD for Kp, 10~ 3 LHD 
for K 1 , and 3 • 10~ 4 LHD for Kg), and, with the excep- 
tion of Kp, practically vanish at LHD. The density de- 
pendence of the yields can be explained qualitatively in a 
way similar to the discussion of the 2s population above. 
At 10~ 7 LHD the cascade is almost purely radiative and 
proceeds at low n through the circular states, so that al- 
most all fi~p atoms go through the radiative 2p — > Is 
transition whereas the transitions np — > Is with n > 2 are 
much weaker. With increasing density the Stark mixing 
becomes more efficient and leads to higher populations of 
the np states with n > 2. This explains the increase in the 
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relative K>p yields and consequently also the decrease in 
K a . The decrease in the K>p yields for higher densities is 
due to the collisional deexcitation n — » n — 1 which begin 
to dominate the radiative transitions. 

At high density, a significant fraction of the fi~ p atoms 
slows down at the very end of the cascade before reach- 
ing the ground state and populate the metastable 2s state 
with high probability (i?2s ~ 0.65 in liquid hydrogen) as 
discussed above. This makes the ratio of K a and Kp yields 
in liquid hydrogen very sensitive to a quenching mecha- 
nism of the 2s state. As demonstrated in Figure [?|, the 
X-ray data favor non-radiative quenching (n4j) m agree- 
ment with the recent experiment on direct search for the 
metastable 2s state | fL2"[ , It would be very important to 
check this result directly by measuring the absolute X- 
ray yield in liquid hydrogen: the non-radiative quench- 
ing leads to a significant suppression of the X-ray yield: 
Y(K tot ) a 0.35. 




Fig. 8. The calculated cascade time and the partial cascade 
time in fi~p vs. density. The partial cascade times are shown 
for different values of the parameter S. The experimental data 
are from |12|. 



3.1.3 Cascade time 

The cascade time in muonic hydrogen at low pressures 
(0.25-16 mbar) was measured recently in the (n~p)i s dif- 
fusion experiment at PSI Jl2] |. The experimentally mea- 
sured value is the so-called partial cascade time that is 
defined as the time between formation of the fi~p and 
the last "significant" collision during the cascade, and is, 
therefore, smaller than the cascade time usually defined 
as the time between formation and arrival in the ground 
state. The exact definition of the measured partial cascade 
time is contained in a non-trivial way in the data analy- 
sis program of [Q; we can approximate it by defining a 
collision to be "significant" if the relative change in pT p 
kinetic energy is larger than 5 which we vary in the range 
0.1 — 0.3. The calculated cascade time and partial cascade 
time are shown in Figure || in comparison with the exper- 
imental data Jl2| . Whereas the cascade time grows as the 
pressure decreases and eventually reaches the purely ra- 
diative cascade time, the partial cascade time reaches its 
maximum at 2 • 10~ 7 LHD and vanishes at zero pressure 
where there are no collisions. The calculated partial cas- 
cade time depends only weakly on S, and our results are 
in good agreement with the experiment. 

The n dependence of the partial cascade time at 0.25 
and 4 mbar is shown in Figure 0. At 4 mbar, the partial 
cascade time depends only weakly on the initial n and 
is in perfect agreement with the experimental result for 
n > 10. The calculated partial cascade time at 0.25 mbar 
has a moderate dependence on n and agrees with the ex- 
periment for n = 13 — 15. Though this seems to support 
the simple picture of the p~p atoms being formed with 
n ~ 14, the models that predict lower values for the ini- 
tial n can also be consistent with the experiment because 
the n distributions usually have tails that extends up to 
high values. For example, using the the n distribution for 
molecular target from we obtain the partial cascade 
time T casc = 60 ns at 0.25 mbar in good agreement with 
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Fig. 9. The partial cascade time as a function of the initial n 
at 0.25 mbar and 4 mbar. The average initial energy is 0.5 eV 
and the parameter 8 — 0.1. The experimental data are from 
0. 



3.1.4 Kinetic energy distribution 

The kinetic energy distribution of muonic hydrogen atoms 
on arrival in the Is state has been obtained by analyzing 
fi~p diffusion times in hydrogen gas |g,0,^,g£]. The in- 
tegrated kinetic energy distribution^ of muonic hydrogen 
on arrival in the Is state 



W(T) 



w(T')dT' 



(15) 



was determined recently at PSI |T^, ^9[. The results ob- 
tained in the pressure range 0.06 — 16 mbar offer a unique 
possibility to study the formation of the p~p and the 



the measured r ri 



74 ± 15 ns \r. 



4 The reason for using the integrated instead of the differ- 
ential kinetic energy distribution, w(T), is that neighboring 
energy bins are highly correlated in the experimental analysis. 
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Fig. 10. The integrated energy distribution W(T) of the fi p 
atom at the end of the cascade for initial conditions: rii — 14 
and T = 0.5 eV. The data are from 111. 



initial stages of the cascade. Figure [H] shows the inte- 
grated kinetic energy distribution calculated in the cas- 
cade model and the experimental data from . The data 
at 0.06 mbar are sensitive mainly to the initial energy dis- 
tribution. With the density increasing, the role of the colli- 
sional processes grows, and the final distribution becomes 
more energetic due to acceleration at the beginning of the 
cascade. The results of the cascade calculations are in a 
fair agreement with the data for the initial conditions spec- 
ified in Section 2.3. The Coulomb transitions with An > 1 



Fig. 11. The integrated energy distribution W(T) of the /i~p 
atom at the end of the cascade for initial conditions: rn = 14 
and To = 0.5 eV. The classical-trajectory calculations at n > 7 
were done for atomic target. The data are from | Q . 



were found to be essential for explaining the observed in- 
crease in kinetic energy with increasing pressure in agree- 
ment with the earlier analysis |ll[] . Similar cascade calcu- 
lations using the trajectories of u~ p collisions with atomic 
hydrogen, where the collisions with An = 1 dominate the 
collisional deexcitation, predict a much weaker increase in 
kinetic energy as shown in Figure [H]. 

The cascade calculations predict, independent of the 
initial n, significantly smaller fractions of (u~p)i s atoms 
with high kinetic energies (T > 8 eV) at pressures above 
4 mbar than the experimental data. 

Figure [j^ shows the density dependence of 14^(1 eV) 
and W(8 eV). The calculated sub-1 eV fraction, W(l eV), 
decreases from 80% at 0.06 mbar to 30% at 16 mbar in 
good agreement with the experimental data p2| . The fi~p 
atoms with energies above 8 eV are produced during the 
cascade in Coulomb transitions. The present calculations 
predict an increase in W(T > 8 eV) with the pressure 
from 1% at 0.06 mbar to 17% at 16 mbar. The predicted 
curves agree with the experimental results in the range 
0.06 — 1 mbar whereas the measured high energy fraction 
is substantially larger than calculated one for pressures 
above 4 mbar: for example W{T > 8 eV) = (22.5 ± 0.9)% 
at 16 mbar JT2| . The measured increase in the high energy 
fraction of (fi~p)i s atoms cannot be reproduced in the 
current cascade model for any initial distribution. 

Figures [l3|, |l4|, and 15 show the density dependence of 
the median kinetic energy on arrival in the ground state 




Fig. 12. The density dependence of the integrated energy 
distributions W(l eV) and W(8 eV) of the ii~p atom for initial 
conditions: rii = 14 and To = 0.5 eV. The data are from 621. 



for different initial conditions. With the density increas- 
ing the calculated median energies grow and reach 2.5 eV 
around 10~ 5 LHD. Above 10 -5 LHD, the median energies 
remain nearly constant up to 2 ■ 10~ 4 LHD where they 
start to grow again. The shape of the curves can be ex- 
plained as follows. The increase in the median energies at 
low densities is produced in the classical domain (n > 7) 
where the acceleration due to Coulomb transitions with 
An = 1 — 5 is more efficient than the slowing down. The 
plateau is due to the dominance of Auger deexcitation 
at medium n which become important around the den- 
sity 10 -5 LHD. The increase in the median energy in the 
uppermost part of the shown density range is produced 
by the Coulomb transitions at n < 8. For densities larger 
than 10~ 4 LHD, the initial kinetic energy is almost forgot- 
ten because of the many collisions during the cascade, and 
the kinetic energy distribution at the end of the cascade is 
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Fig. 13. The density dependence of the median kinetic energy 
of the fi~p at the end of the cascade for different initial average 
kinetic energies and m = 12. The data are from 
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Fig. 14. The density dependence of the median kinetic energy 
of the fj,~p at the end of the cascade for different initial average 
kinetic energies and rn = 14. The result for atomic target is 
shown for T = 0.5 eV. The data are from @G3. 
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Fig. 15. The density dependence of the median kinetic energy 
of the fi~p at the end of the cascade for different initial average 
kinetic energies and rii = 16. The data are from jl^,Q. 



Observable 


ESCM 


Exp. 


T mcd ian(l mbar) 


1.57 ±0.03 eV 


1.55 ±0.12 eV 


T mc dian(16 mbar) 


2.50 ±0.14 eV 


2.63 ± 0.17 eV 


W(16 eV) at 16 mbar 


0.958 ± 0.004 


0.894 ± 0.006 



Table 2. Median kinetic energy and integrated kinetic energy 
in muonic hydrogen calculated in the ESCM with statistical 
errors. The experimental results are from |J12[|44|. 



does not change the conclusions reached above: the ESCM 
gives a good description of the experimental data for mod- 
erate energies but is unable to explain the observed high 
energy fraction. 



3.2 Antiprotonic hydrogen 

The present study of the atomic cascade in antiprotonic 
hydrogen was focused on the experimental data obtained 
at low density &|f|. 



determined by the competition between acceleration and 
deceleration during the cascade. 

The measured density dependence of the median ki- 
netic energy is well reproduced with the initial conditions: 
T = 0.5 eV andn, = 14 or T = 0.25 eV and m = 16. For 
rii = 12 (Figure |l3|) the initial acceleration appears to be 
weaker and the agreement with the experimental result is 
worse. The molecular structure of the target is essential 
for explaining the data. Calculations with the CMC tra- 
jectories for atomic hydrogen leads to a too weak increase 
in the median energy compared to the experiment. 

We estimate the statistical errors due to the finite num- 
ber of CMC trajectories by dividing the 21 complete sets 
of trajectories into 7 subsets and calculate the observables 
using each subset. The results for the median kinetic en- 
ergy and the integrated kinetic energy shown in Table ^| 



3.2.1 X-ray yields 

The calculated X-ray yields are in a good agreement with 
the data [^6|— 149|] as shown in Figures H and [L7|. In addi- 
tion to the data on absolute X-ray yields, there are more 
precise measurements of relative yields, see Table |[ The 
theoretical results are in fair agreement with the data as 
all the calculated ratios deviate less than 50% from the 
experimental results. This is a significant improvement in 
comparison with the calculations in the Borie-Leon model 
reported in |^| : the two standard tuning parameters fcstk 
and T were fixed by fitting the absolute L a yields but the 
predictions of the model for ratios involving L$ were a fac- 
tor of 3 — 5 larger than the experimental data. We stress 
that our results were obtained without employing any cas- 
cade tuning parameters. Figure [l^ shows an example of 
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Fig. 16. The density dependence of the absolute L X-ray Fig. 18. The K line profile in antiprotonic hydrogen at 
yields_in antiprotonic hydrogen. The experimental data are 10 -3 LHD. 
from 
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Fig. 17. The density dependence of the absolute K X-ray 
yields in antiprotonic hydrogen. The experimental data are 
from 



Fig. 19. The density dependence of the cascade time in an- 
tiprotonic hydrogen. The experimental data are from p5[. 



Rad. 



the X-ray profile of the K lines at 10~ 3 LHD. Because of 
the hadronic broadening of the Is state, only the K a line 
is clearly separable from the rest of the K lines of which 
the 9p — » Is transition has the highest intensity. The rela- 
tively large contribution to the K yield from high n states 
is due to the strong absorption in the excited states which 
prevents most of the pp atoms from reaching the low n 
states. 
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3.2.2 Cascade time 
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The density dependence of the cascade time is shown in 
Figure [l9[ Our result is in good agreement with the data 
p5fl except for the data point at 150 mbar where our model 
predicts significantly shorter cascade time. The agreement 
for pressures up to 10 mbar can be considered as an impor- 
tant confirmation of the dominance of the Coulomb transi- 
tions with large changes in the principal quantum number 



Fig. 20. The dependence of the cascade time in pp at 3.4 and 
9.8 mbar on the initial state n. The initial average kinetic en- 
ergy is 0.5 eV. The cascade time for a purely radiative cascade 
is shown for comparison. The experimental data are from [M . 
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Observable 


Pressure 


ESCM 


Exp. 


Reference 


Y{K a )/Y(L a ) 


20 mbar 


0.0179 


0.0176 ±0.0016 


1 


51 
51 

48 
48 
48 
48 
48 
48 
48 




Y{L tot )/Y(L a ) 


20 mbar 


1.56 


1.45 ±0.05 


I 




Y(L a )/Y(L p ) 


30 mbar 


5.0 


7.1 ±0.7 


I 




Y(L a )/Y(L^) 


30 mbar 


13.2 


20.6 ± 2.0 


I 




Y{L a )/Y(L s ) 


30 mbar 


22.6 


23.9 ± 2.6 


I 




Y(L a )/Y(L tot ) 


30 mbar 


0.62 


0.75 ±0.08 


I 




Y(L P )/Y(LJ 


30 mbar 


2.7 


2.9 ±0.2 


I 




Y(L I3 )/Y(L S ) 


30 mbar 


4.6 


3.4 ±0.2 


I 




Y(L~ f )/Y(L s ) 


30 mbar 


1.7 


1.2 ±0.1 







Table 3. Relative X-ray yields in pp. 




N/LHD 

Fig. 21. The density dependence of the median energy at 
the instant of nuclear absorption in antiprotonic hydrogen for 
initial n = 20, 25, 30. 



n in the upper part of the atomic cascade. The earlier 
calculations || predicted much higher cascade times in 
variance with the experimental data. 

The dependence of the cascade time at 3.4 and 9.8 mbar 
on the initial n is shown in Figure ^0. There is good 
agreement with the experimental result [l5| in the range 
n = 25 — 35. Calculations also show that the calculated 
cascade times depend only weakly on the initial kinetic 
energy. 
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3.2.3 Kinetic energy distributions 



Fig. 22. The calculated integrated energy distribution of pp 
at the instant of the 3d — > 2p transition at 22 mbar. 



3.2.4 Doppler broadening of the L lines 

The acceleration processes during the atomic cascade lead 
to a sizeable Doppler broadening of the L lines. Figure [2^ 
shows the integrated energy distribution of the pp atoms 
at the instant of the 3d — > 2p transition in a gaseous target 
of 22 mbar, the median energy being about 6.5 eV. The 
corresponding Doppler broadening of the 3d — > 2p line 
profile is shown in Figure [2^; the effect gives a significant 
contribution to the width and must, therefore, be taken 
into account in the analysis of the X-ray spectra. In the 
data analysis in |Q it was assumed that the pp atoms 
were thermalized due to elastic and Auger collisions. 



In the upper part of the cascade most of the pp atoms 
are accelerated to energies of several eV in the Coulomb 
transitions with An > 1. Figure pl| shows the density de- 
pendence of the median energy at the instant of nuclear 
absorption calculated for initial n = 20, 25, 30. The shape 
of the curves are similar to those of muonic hydrogen but 
the energy is higher. The median energy increases rapidly 
with increasing density and reaches 6 — 7 eV at 10 -5 LHD. 
At higher densities the dominance of the Auger transi- 
tions at intermediate n prevents further increase until the 
Coulomb deexcitation at low n becomes important. 



4 Conclusions 

Using the new results for the collisional processes 
the atomic cascades in muonic and antiprotonic hydrogen 
have been have been studied in detail. The cascade cal- 
culations have been done in the extended standard cas- 
cade model which describes the evolution in the quantum 
numbers n and I and the kinetic energy from the very be- 
ginning of the cascade. By taking deceleration and accel- 
eration during the cascade into account, observables can 
be predicted reliably without the need for any cascade 
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Fig. 23. The natural L a line profile in pp (dashed line) in 
comparison with the Doppler broadened line profile at 22 mbar 
(solid line) . For the electromagnetic hyperfine splitting we use 
the result of Boucard and Indelicate [pOj. The hadronic split- 
ting and broadening are taken from the DR1 model of I 



the high energy fraction of the kinetic energy distribution 
in muonic hydrogen deserves special attention in future 
studies: the ESCM predicts a much smaller fraction of n~p 
atoms with energies above 16 eV at 4 and 16 mbar than the 
experimental result |i~2] ] . This disagreement indicates that 
some processes beyond the ESCM are important. Can- 
didates under discussion are the hybrid Auger-Coulomb 
transitions [|53f and the resonant formation of molecular 
states f§. 

Future development of the ESCM will consist of calcu- 
lations of the above mentioned processes and their inclu- 
sion in the cascade model. At the same time the cascade 
calculations can be easily extended to cover exotic deu- 
terium atoms and H-D mixtures. Further critical tests of 
the current theory of the atomic cascade will greatly ben- 
efit from the results of the experiments on the precision 
spectroscopy of pionic and muonic hydrogen, presently in 
progress at PSI. 
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tuning parameters. The extended standard cascade model 
successfully describes a number of experimental results: 
X-ray yields, cascade times, and kinetic energy distribu- 
tions. The results of the cascade calculations for muonic 
and antiprotonic hydrogen can be summarized as follows: 

— Muonic hydrogen. The calculated X-ray yields in muonic 
hydrogen are in good agreement with the data. The 
measured ratio of the K a and Kp yields in liquid hy- 
drogen can be reproduced only if the non-radiative 
quenching of the metastable 2s state is dominant. The 
calculated increase in the ([i~p)i s median kinetic en- 
ergy with the density is in agreement with the data |jl2f 
for the initial conditions corresponding to kinetic en- 
ergies around 0.5 eV and principal quantum number 
n 14. The molecular structure of the target is essen- 
tial for explaining the data on the density dependence 
of the kinetic energy distribution at the end of the cas- 
cade. 

— Antiprotonic hydrogen. The calculated K and L X-ray 
yields in antiprotonic hydrogen are in good agreement 
with the data. The calculated cascade time is in good 
agreement with the data. For the first time, a good 
description of the data is obtained without using any 
phenomenological tuning parameters. The fast colli- 
sional deexcitation at high n calculated in the classical- 
trajectory model was found to be very important for 
explaining the short cascade times measured by the 
OBELIX collaboration |^] . The acceleration of the pp 
during the initial stages of the cascade predicted by our 
cascade model results in a significant Doppler broad- 
ening of the L lines. This must be taken into account 
when hadronic widths and shifts are obtained from the 
measured line profiles. 

Some experimental data cannot be explained in the 
present version of the ESCM even with taking the uncer- 
tainties in the cross-sections into account. The problem of 
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